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The magnetic flux dependence of order parameter and supercurrent is studied in a hollow d-wave 
superconducting cylinder. It is shown that the existence of line nodal quasiparticles in a pure d x z_ y 2 
pairing state gives rise to an hc/e periodicity in the order parameter and a first-order quantum phase 
transition for a large system size. We demonstrate that the flux periodicity in the supercurrent is 
sensitive to the detailed electronic band structure and electron filling factor. In particular, we find 
that, in cooperation with the increase of the cylinder circumference, the hc/2e periodicity can be 
restored significantly in the supercurrent by avoiding the particle-hole symmetry point. A similar 
study of a d x 2_ y 2 +id xy pairing state verifies the peculiarity of unconventional superconductors with 
nodal structure. 
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A fundamental property of all known superconductors 
is the formation of Cooper pairs [l[ in the superconduct- 
ing state. A far-reaching implication of this fact is the 
quantization of magnetic flux in units of hc/2e in mul- 
tiply connected superconducting geometries. The hc/2e 
flux quantization has been used as a proof of electron 
pairing nature of both conventional @, y, 0, [j| and high- 
temperature Q superconductors in the superconducting 
state. Other related phenomena include the quantum os- 
cillation in the transition temperature Q and magnetic 
vortices each carrying an hc/2e flux quantum [1, 

Quantum mechanically, there is no fundamental rea- 
son why the minimal flux periodicity must be hc/2e in 
a superconductor. The gauge invariance can only guar- 
antee a fundamental period of <& = hc/e [3, EH- Only 
when all Cooper pairs move in the same group veloc- 
ity, can a substantial $o/2 periodicity be obtained. The 
$o-periodicity has been in mesoscopic conventional su- 
perconducting rings [111, LL2|, [l3|, ll4[ due to the level dis- 
creteness and Landau depairing effect. More surprisingly, 
recent studies have shown a severe breakin g o f hc/2e- 
periodicity in a d-wave superconducting loop |H,[H,[l3|; 
as a result of the Cooper-pair angular momentum se- 
lection for the existence of Doppler-shifted zero-energy 
states. 

In this Letter, by a systematic study of the supercur- 
rent in a hollow <i-wave superconducting cylinder, we 
show an hc/e magnetic periodicity in the d-wave order 
parameter and demonstrate that the flux periodicity in 
the supercurrent is sensitive to the detailed electronic 
band structure and electron filling factor. In particu- 
lar, we find that the breaking of hc/2e periodicity in the 
case of d-wave pairing is closely related to the particle- 
hole symmetry in the normal state band structure, which 
gives rise to the van Hove singularity. When the particle- 
hole symmetry point is avoided, the hc/2e periodicity can 
be restored almost completely in the supercurrent. 

To be specific, we consider a hollow d-wave super- 
conducting cylinder, as schematically shown in Fig. [1] 




FIG. 1: (Color online) Schematic drawing of a hollow rf-wave 
superconducting cyliner. A magnetic flux threads the cylinder 
parallel to its axis. 



Experimentally, the cylinder can be formed by a high- 
temperature superconductor film with its normal perpen- 
dicular to the Cu02 plane. A magnetic flux <& threads 
the cylinder parallel to its axis, and also to the crystal 
6-axis of the Cu02 plane. Due to the weak interlayer 
coupling, the problem can be reduced to one on an es- 
sentially two-dimensional system. We define the x- and 
y-axis to be perpendicular and parallel to the flux di- 
rection, respectively. This set up can also avoid the nu- 
clcation of Abrikosov vortices, which will complicate the 
analysis. By choosing a gauge, where the vector poten- 
tial does not appear explicitly in the Hamiltonian, the 
Bogoliubov-de Gennes (BdG) equations can be written 
as 1181: 
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Here are the eigenfunctions corresponding to 

eigenvalues E n . The single particle Hamiltonian H.^ = 
—tij — fiSij with tij and /i being the hopping integral 
and chemical potential. We consider the nearest neigh- 
bor, t, and next-nearest neighbor, t' , hopping integral. 



The bond order parameter for (i-wave pairing is deter- 
mined self-consistently as Ay = (V x 2_ y 2 /4) J2 n l u f v j * + 
Uj J v™*]t&nh(E n /2kBT) with V x 2_ y 2 being the pairing 
strength in the d x 2_ y 2 channel. Notice that the quasi- 
particle excitation energy is measured with respect to 
the Fermi energy. 

If the electron wave vector is k = (k x ,k y ) and that 
for the collective drift motion (supcrfluid motion) of the 
paired electrons is q = (q x ,q y ), the initial k and k 
pairing is adjusted to pair the states k + q and k + q. 
The solution to the BdG equations is then found as: 
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Here, corresponding to E^ = Zk, q ± £? k " ' the electron 
and hole components of the quasiparticle amplitude are 
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FIG. 2: (Color online) Flux dependence of the persistent cur- 
rent for a normal state metallic cylinder of size Nl = 40 2 
(solid line) and Nl = 80 2 (dashed line), with n e = 1.0 and 
t' — 0. The current is measured in units of To — ei. The other 
parameter values are defined in the main text. 
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The quantities Qk.c 



[£k+q + 6c-q]/2, ^k.q — [6c+q - 

[Qk, q + A k] 1/2 - In thc tig 11 *- 

binding approximation, up to the next-nearest neighbor, 
the conduction electrons have the normal state disper- 
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The d-wave superconducting gap dispersion is given by 



Ak = 2A,i<^k, where 0k = cosfcj 
determined self-consistently: 



cosfcj,, and A^q) is 
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where A^ = N x x TVy and the quantities N x and A^ rep- 
resent the circumference and length of the hollow cylin- 
der. Rigorously, the bond order parameters along the x 
and y directions do not follow the relation A x = —A y 
in thc presence of magnetic flux. Here we have imposed 
the restriction of A x = — A y to enforce a rigorous d- 
wave symmetry. Notice that A^ is now a function of 
From the boundary condition given by Eq. ([2]), 
one can find [lH that the components of wave vectors 
k x = 2ir(n x -m/2)/N Xl q x = 2<K{$/§ +m/2)/N x , while 



k y = 2nn y /Ny, and q y 



0, where 



l x(y) 



and m arc inte- 



gers. In particular, m is determined by minimizing \q x \ 
for a given value of magnetic flux <!>. The electron fill- 
ing factor and the single nearest-neighbor bond current 
flowing around the cylinder are computed via, 
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respectively A factor of 2 has been included to account 
for the spin degeneracy 

In thc numerical calculations, we take k B = t = 1. 
Throughout the work, the energy is measured in units of 
t unless specified otherwise. The temperature is fixed at 
T = 0.01 and the d-wave channel pairing interaction is 
chosen to be V x 2_ y 2 = 4. Both the hopping parameter t' 
and thc electron filling n e will be changed. For a given n e , 
the chemical potential should be adjusted and, therefore, 
will be a function of $/$o- 

In Fig. [21 we show the flux dependence of the persistent 
current in a normal state metallic cylinder, where there is 
no existence of superconducting order parameter. These 
results are known from the study of persistent current in 
normal state mesoscopic rings hjl in the presence of an 
Aharonov-Bohm flux [2(J. The main point is that the 
persistent current in a normal state ring has a period- 
icity of <i>o- We present them here to demonstrate that 
our theoretical formulae designed for the superconduct- 
ing state can reduce to describe the normal state, and to 
provide a starting point for our discussion on the case of 
d-wave superconducting state below. 

Figure [3] shows the flux dependence of thc e?-wavc order 
parameter and persistent current for various electron fill- 
ing and system size but with t' = 0. Due to the existence 
of nodal quasiparticle in a d x 2_ y 2 pairing symmetry, the 
order parameter oscillates in the magnetic flux with a pe- 
riod of <&o, even for larger cylinders, although the oscil- 
lation amplitude is even smaller. The oscillation pattern 
is strongly sensitive to the system size. In particular, it 
is seen that the system experiences a first-order quantum 
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FIG. 3: (Color online) Flux dependence of the d-wave order 
parameter (al-a2) and persistent current (bl-b2) for a rf-wave 
superconducting cylinder of zie Nl = 40 2 (al-bl) and Nl = 
80 2 (a2-b2), with n e = 1.0 (solid line) and n e = 0.8 (dashed 
line). In (al-a2), the relative amplitude of the d-wave right 
parameter, [Ad($) — A<j(<1? = 0)]/A<j($ = 0), is shown. Here 
i! = and the other parameter values are defined in the main 
text. 



phase transition when the flux across 3>/$o = (2n — l)/4 
with n an integer, by exhibiting a discontinuous jump. 
The jump amplitude also decreases with increased system 
size, indicating the first order quantum phase transition 
will be changed into a second order one in the thermo- 
dynamic limit. The $o periodicity in the pairing order 
parameter renders that the supercurrent does not develop 
fully the $o/2 periodicity, in a mathematically rigorous 
sense, regardless of the system size and electron filling. 
The reason lies in the fact that for a pure d x 2_ y 2-wa,ve 
superconductor, low energy quasiparticlc states can be 
populated by the Doppler shift rapidly along the nodal 
direction, for which the orientation dependent coherence 
length is divergent. As a hallmark of superconducting 
state, one can see clearly that flux-induced current has a 
different sign from that of the normal state in the region 
approximately from $/$ = n — 1/4 to <&/$o = n + 1/4 
(compare Fig. |3fbl-b2) with Fig. 0). At the half fill- 
ing, the particle hole symmetry holds, which makes a 
large level spacing at the fermi energy when the circum- 
ference of the cylinder is small. Therefore, the super- 
current exhibits an activation-like behavior for the mag- 
netic flux close to $/$o = n, an d is very different from 
that at $/$o = n/2. It explains why the 4>o period- 
icity of supercurrent is pronounced in this specific case. 
Furthermore, it is not difficult to anticipate that, as a 
canonical mesoscopic effect, the breaking of hc/2e pe- 
riodicity will be even more stronger for a smaller zero- 
field e?-wavc pair potential, and therefore a larger momen- 
tum averaged superconducting length £ = (hv^/nA^ps, 
where is the quasiparticlc velocity. We notice that, 
in the region <f>/$o G [n - 1/4, n + 1/4], the flux de- 



FIG. 4: (Color online) The Same as Fig.[3]but with t' = -0.2. 



pendence of the supercurrent in a d-wave superconduct- 
ing cylinder is different from that in a square d-wavc 
superconducting loop [l5|, where a zig-zag feature was 
obtained. We argue that in the geometry considered in 
Ref. [TBI, the elastic scattering from hard-wall boundaries 
of a mesoscale system populates a significant number of 
the lower energy states, which plays an important role 
in the flux-dependent current. Naturally, the increase of 
the cylinder circumference is one way to reduce the acti- 
vation behavior, and therefore reducing the $/$o com- 
ponent in the Fourier spectrum of supercurrent. Alter- 
natively, when the electronic filling factor is tuned away 
from the half filling, the particle-hole symmetry is bro- 
ken and the Fermi surface becomes more isotropic. The 
activation-like behavior in the supercurrent is replaced by 
a more linear-like behavior, similar to the behavior ex- 
hibiting at $/$o = n/2. In addition, the current peaks 
at $/$o = (2n — l)/4 becomes more symmetrized about 
1 = axis, tending to restore the more of $0/2 per- 
odicity in supercurrent. Figure U shows the flux depen- 
dence of the d-wave order parameter and persistent cur- 
rent for various electron filling and system size but with 
t' = —0.2. When a finite next nearest-neighbor hopping 
integral is introduced, the particle-hole symmetry is bro- 
ken at the outset for the hole doped (i.e., n e < 1) system, 
where the chemical potential is not zero. In this case, the 
level repulsion at the Fermi energy is weakened for flux 
close to $/$o = n even for a small cylinder circumfer- 
ence, and the activation behavior in the current does not 
show up. Consequently, the $0 component in the current 
spectrum is dramatically decreased, which makes the to- 
tal current looks to have $o/2 periodicity completely. We 
note that when the zero-field averaged superconducting 
coherence is not so small in comparison to the cylinder 
circumference, the restoration is always incomplete, due 
to the mesoscopic effect. 

To understand better the magnetic flux periodicity in 
unconventional superconductors, we turn to consider a 
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flux dependence of the order parameter similar to that 
shown in Fig. 2 in Ref. 11 or Fig. 6 in Ref. [ID, which 
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FIG. 5: (Color online) Flux dependence of the ci-wave or- 
der parameter (al-a2) and persistent current (bl-b2) for a 
d x 2_ y 2 +id xy -W8<ve superconducting cylinder of zie Nl = 40 2 
(al-bl) and N L = 80 2 (a2-b2), with n e = 1.0. In (al-a2), the 
relative amplitude of the respective d x 2_ y 2 (solid line) and 
d xy (dashed line) components are plotted. Here t' = and 
the other parameter values are defined in the main text. 



d x 2_ y 2 +id xy pairing state by taking the pairing strength 
in the id xy channel as V xy = 3.0. Now the quasiparticle 
excitions are gapfull. In Fig. we show the flux depen- 
dence of the d-wave order parameter and supercurrent 
in a hollow d x 2_ y 2 + id xy -wave superconducting cylin- 
der. Noticeably, even in the presence of the particle-hole 
symmetry and for the same system size as the case of 
a pure d x 2_ y 2 pairing state, both the d x 2_ y 2 and d xy 
components of order parameter have the periodicity of 
$o/2 = hc/2e. In particular, one case see that the evolu- 
tion of two components of order parameter is continuous 
when $/$o crosses (2n — l)/4, indicating that the flux- 
induced first order quantum phase transition is unique 
to a cylinder formed by unconventional superconductors 
with nodal quasiparticles. In the present case, the mag- 
netic hc/2e periodicity in the current is complete. The 
hc/2e periodicity is set in as long as the cylinder circum- 
ference is much larger than the superconducting coher- 
ence length. We point out (but do not show) that the 
flux dependence of order parameter and supercurrent in 
a hollow s-wave superconducting cylinder exhibit similar 
behavior to the case of d x 2_ y 2 + id xy pairing state. 

One remark is in order: Our calculations have shown 
that the magnitude of the order parameter can be en- 
hanced in the presence of magnetic flux throughout the 
whole period of 4> (see e.g., Fig. [3fa2) and Fig. [5fal)- 
(a2)). This is the characteristic of all superconduc- 
tors (including s-wave case) with short coherence length, 
which is much smaller than the cylinder circumference. 
We have calculated the case of a pure d x 2_ y 2-wave su- 
perconductor and a conventional s-wave superconductor 
both with a pairing interaction equal to 1 , and found the 



again is a mesoscopic effect. 

In conclusion, we have studied the flux dependence of 
order parameter and supercurrent in a hollow d-wave su- 
perconducting cylinder. For a pure d x 2_ y 2 pairing state, 
we find a hc/e periodicity of order parameter due to the 
existence of nodal quasiparticle states, and an associated 
quantum phase transition. When the particle hole sym- 
metry holds in the normal state band structure, there is 
a noticeable component of hc/e in the supercurrent spec- 
trum when the cylinder circumference in the mesoscopic 
regime. However, in addtion to the increase of system 
size, this component can be suppressed more effectively 
by avoiding the particle hole symmetry point through the 
tuning of electron filling and band structure. By study- 
ing the case of a d x 2_ y 2 + id xy pairing state, where the 
quasiparticle excitations are gapfull, we verify that the 
peculiar hc/e magnetic flux periodicity only happens to 
unconventional superconductors with nodal structure. 
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